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We consider the propagation of light from an isolated source that also emits gravitational waves. 

The eikonal approach is employed to determine the transfer of energy from the gravitational to 
electromagnetic radiation. A mechanism is found in which a photon “surfs” on the gravitational 
wave. For black hole events, a signihcant upshift of photon energy can occur according to a power- 
law buildup over the radial distance. This surprising effect may be responsible for some of the 
unexplained high energy phenomena in the cosmos involving gamma rays or other astro-particles. 
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The ray description of light has long offered a powerful 
theoretical tool in probing spacetime structure and its ef¬ 
fect on matter kinetics. It is the very approach adopted 
by Einstein, together with the constancy of the speed 
of light, in arriving at the celebrated special relativity a 
century ago. The subsequent advent of general relativity 
requires far more sophisticated analysis of null geodesics 
and light cones to unravel local and global spacetime 
properties that are often counterintuitive. These include 
the existence of black hole and cosmological horizons as 
well as gravitational lensing, focusing and red/blue shifts. 

The anticipated direct detection of gravitational waves 
(GWs) 0 and ongoi^ intense research into GW inter¬ 
action with matter 00 ; 0 and astro-particle physics 
0 have stimulated a great surge of interest in novel as¬ 
tronomical phenomena that involve both high energy jets 
and gravitational radiation 00 For events as violent as 
the core-collapse of stars and merger of black holes, one 
of the largest portions of energy released is through grav¬ 
itational radiation. This leads to a number of important 
issues. For instance, does there exist any mechanism by 
which some of the GW energy may be transferred to pho¬ 
tons and other particles? If such a mechanism exists, how 
relevant is it for energetic events like gamma ray bursts 
(GRBs) 000? Furthermore, to what extent can such 
a mechanism be responsible for the unexplained sources 
of energy for the observed ultra-high energy cosmic rays 

0 ? 

In this letter, we address these issues by exploring a 
hitherto undiscovered process in which light particles and 
GWs from an unstable compact source may interact in a 
resonant manner, causing significant energy gain by the 
particles. In particular, we shall focus on the GW in¬ 
teraction with photons with a view to further dedicated 
treatments for other particles such as neutrinos. The ray 
description of light propagation will be adopted, as ap¬ 
propriate for photons (at 1 eV or larger) having a wave¬ 
length much shorter than that of the GWs (of a few km 
or larger) under consideration. 


In some sense, the effect of photon energy upshift 
W GWs is in close analogy with that by plasma waves 
[SUlflj. In the latter case, as a medium is involved, pho¬ 
ton energy upshift is accompanied by the acceleration 
of the photon group velocity, the effect has been coined 
as “photon acceleration”. Some pioneering theoretical 
works have been done to replace plasma waves by ide¬ 
alized plane GWs 00 , where it has been found that 
the photon energy may increase for a period before it 
falls. Recent development of the “collapsar” model of 
GRBs strongly suggests that gamma ray bursters may 
emit significant GWs 000. However, to date no at¬ 
tempts have been reported to couple the predicted GWs 
with the gamma rays. In view of this, the present work 
will consider the gravitational radiation and explore the 
consequent interaction with co-moving photons. 

Below, the geometric units of c = G = 1 are adopted 
unless otherwise indicated. The notation for differential 
geometry follows mainly that of 0 - In particular, the 
Greek and Latin coordinate indices range over 0-3 and 1- 
3 respectively. The essential picture to envisage consists 
of an isolated source of GWs and photons. The former 
perturb the asymptotically fiat spacetime such that, in 
coordinates {x^) = = {t,x,y,z), the spacetime 

metric assumes the form 

fjpi,' ~ “f (^) 

where 77 ^,^ = diag(— 1 , 1 , 1 , 1 ) and > 0 at infinity. 

This paradigm is appropriate for astrophysical conditions 
where the expansion of the Universe is negligible, as will 
be justified towards the end of this letter. The total 
energy of the photons are taken as insufficient to back- 
react on the GW metric. The radiation from the source 
will be detected by a remote observer. Both the source 
and observer are assumed to be stationary with respect 
to the Minkowski metric and follow the asymptotic 
“cosmological time” flow given by d/dt. 

Note that although there exist a large class of 
Robinson-Trautman solutions for spherical GWs, the as- 
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sociated singularity structures unfortunately render these 
solutions unsuitable to described gravitational radiation 
from an isolated source In astrophysics, GWs from 
such a source have very small amplitude in the far field, 
where the GW interaction with photons under consid¬ 
eration take place. Therefore, the lowest order post- 
Newtonian approximation to general relativity in terms 
of linearized gravity provides an adequate description. 

In treating the photons, we shall take into account not 
only their trajectories but also their energy in reference 
to the “cosmological time” t = x^. To this end, consider 
the eikonal function ijj = along the light rays. It 

gives rise to the wave 4-vector as = d'lli/dx^^. For any 
spacetime metric g^’', the eikonal function ip satisfies the 
dispersion relation for photons given by g^'^k^ki, = 0. 
The angular frequency of the photon with respect to t 
is given by w := —fco = —dipjdt. (Here means “is 
defined as”.) This may be compared with the proper 
frequency lOt given by Wr := —drip = uj/^y—goo with 
respect to the proper time r along dt such that dr^ = 
—goodt^. It follows from the dispersion relation that the 
positive photon frequency uj takes the form 


iO = 


ujf + UJ 2 — W2 


where 


oji := 



UJ2 ■ = 


g^°’ka 

-gOO- 


( 2 ) 

( 3 ) 


The photon dynamics is generated from a variational 
principle using the action integral 


S = 


kfj, da:'' = 



u) dt 


( 4 ) 


where the over dot denotes a t-derivative. This action 
takes a canonical form in terms of the phase variables 
(a;“,fca) and Hamiltonian uj = oj(x°‘,ka,t) specified in 
0, yielding the following canonical equations of motion: 


duj l _ 


( 5 ) 


It is important to note that the Hamiltonian to is nec¬ 
essarily time dependent if the spacetime metric is non¬ 
stationary. 

Return to the metric of the form 0 and treat h^i, per- 
turbatively to its linear order henceforth. For our pur¬ 
pose it suffices to restrict to hoo = 0, hence goo = = 

— 1 (to linear order in which makes t proper and 

simplifies subsequent analysis. Substituting this condi¬ 
tion and 0 into 0 we see that 


( 6 ) 


where := k'^ + ky + k^. Therefore the canonical equa¬ 
tions of motion 0^ reduce to 

^ h^^^kakbk^ + (7) 

k k 2k'^ 

k = ^d,h‘^'>kakb-d,h°’^ka. ( 8 ) 

ZK 

As in standard linearized gravity, consider the trace- 
reversed metric perturbation := — 4 with 

the trace term h = 77 '“'satisfying the Lorentz gauge 
condition h '^= 0. The gravitating source has the full 
quadrupole moment 

Qab{t) = Jpx°‘x^d^x (9) 


where p denotes the mass density. We shall evaluate 
the metric perturbation at a large distance L compared 
with the size of the source. Within the level of accuracy 
adopted, h^i, is given by the quadrupole formula 

Kb = ^ Qab{t - L) (10) 

along with Ko = Ka = 0 ^3. To express the metric 
perturbation itself, it is more convenient to introduce the 
trace-reversed quadrupole moment given by 

Qab{t) = J p{x°'x^ --6abK)d^x (11) 


where := For simplicity the radial motion 

inside the source is assumed to be negligible. Hence, 

hab — -j- Qab{t — L) = — Qab{t — L) (12) 

along with hoo = Ka = 0. Based on the above discussion 
we may now investigate the dynamics of photons on the 
GW background. Of particular interest in our analysis is 
the asymptotic behaviour of photons originated from the 
source. It is clear that in the far field the photon trajecto¬ 
ries are essentially radial. Let us take a photon travelling, 
say, in the cc-direction. Thus we have L{t) = x{t). For 
simplicity, introduce K{t) := kx{t). By retaining only 
the leading contributions in 0 and ©, we seek to find 
the asymptotic behaviours of L{t) and K{t) starting from 
the initial position L = Lq and momentum K = Kq at 
t = 0. This leads to the following equations 

L = 1, k = Y dxKx- (13) 

The first equation above yields immediately L = t Lq. 
To solve the second equation, we calculate that 

2 ■ 2 •• 

^x^xx — Qxx^^ 2 

X 

using m- Neglecting the second term in d and sub¬ 
stituting the remaining part into d we obtain 


dK 

dT 



w = /c — T h^^kakb + hP°'k, 


(15) 
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where the relation L = t + Lq has been used to turn K 
into a function of L and the dimensionless quantity 

Q ■■= -Q.A-Lo) (16) 

has been introduced. It is worth stressing that Q is a 
constant of motion that depends on the phase of GW at 
the initial position of the photon. The argument —Lq 
in oni counts for the retarded time required for GW to 
reach x = Lq at t = 0 from the centre of the source 
where x = 0. The asymptotic momentum for the photon 
is therefore obtained by solving m to be 

/ r \ Q 

/../>■„(-) . , 17 ) 

Remarkably, this expression states that, depending on 
the sign of Q, the asymptotic momentum of the photon 
can either increase or decrease monotonically, as the re¬ 
sult of the photon “surfing” the GW at different phases. 

This statement has to be augmented with certain phys¬ 
ical considerations. If the photon’s momentum drops to 
the extent that its wavelength is comparable to that of 
the GW, then it can no longer be described as a particle 
and must be treated as a wave. For Q > 0, as the photon 
gains more momentum, it behaves even more particle¬ 
like. However, several effects may come into play to cap 
the photon energy upshift. First, the back-reaction from 
the photons will become significant if their total energy 
is comparable to that of the GW. Secondly, the GW am¬ 
plitude drops as 1/L and will cease to dominate over 
surrounding spacetime geometry when it becomes close 
enough to the level of stochastic GW background. Above 
all, upon reaching the energy levels of 10 ^^- 10 ^®eV, indi¬ 
vidual photons travelling along the GW will be scattered 
by ambient photons including those from starlight and 
the cosmic microwave background (GMB), by exchang¬ 
ing pairs, and hence lose the status of being free 

from non-gravitational interactions. 

It would appear from the dispersion relation that 
the speed of the photon differs from the speed of light, 
which has been normalized to unity. This is simply be¬ 
cause the speed of the photon is evaluated here with re¬ 
spect to the auxiliary Minkowski metric The speed 
of the photon in vacuo is, of course, the speed of light, 
when evaluated with respect to the physical curved met¬ 
ric However, this could cause a concern that photons 
might de-phase from the GWs that travel at the speed of 
light with respect to 77 ^^. Nonetheless, a closer inspection 
of this dispersion relation reveals the following important 
properties. Any photon moving in the x-direction at a 
GW phase where h^x = 0 has speed equal to unity. At 
this phase if dxhxx > 0 then the energy of this photon 
will continue to increase while co-moving with the GW 
at the same speed. Photons following not far behind ex¬ 
periencing hxx < 0 will catch up while those travelling 
slightly ahead and experiencing hxx > 0 will slow down. 



FIG. 1: This diagram illustrates how the energy upshift of 
photons by an expanding GW is enhanced by “phase lock-in”. 
As shown, consider a portion of the h^x component of the GW 
far off the radiating source along the x-axis. It is divided into 
four phases, separated by dashed vertical lines, according to 
the signs of hxx and dxhxx- Owing to the positive dxhxx of 
phases II and III, the co-moving photons within these phases 
gain energy, as indicated by upwards arrows. In contrast, 
photons within phases I and IV, where dxhxx < 0, suffer from 
a loss of energy, as indicated by downwards arrows. With 
respect to the auxiliary Minkowski metric, photons (dots with 
right-arrows) within phases I and II, where hxx < 0, travel 
faster than the GW whereas photons (dots with left-arrows) 
within phases III and IV, where hxx > 0, travel slower than 
the GW. As such, photons tend to be bunched between phases 
II and III, where hxx = 0 and dxhxx may well be maximum. 
The bunched photons are indicated by the un-arrowed dots. 


In the long run, these photons tend to be bunched around 
the GW phases where hxx = 0 and dxhxx > 0 thereby 
enjoying continuous energy rise. This is illustrated in 
Fig.n The same argument applies for photons radiating 
in any other direction. The result is a modulation of the 
photon distribution on the GW wavelength scale. 

Let us remark on the role of gauge at this point. 
Within the framework of linearized general relativity us¬ 
ing the Lorentz gauge condition (without loss of gen¬ 
erality), it is well known that under any infinitesimal 
coordinate transformation with the generating functions 
= ,c^(x'") satisfying = 0 , the metric perturbation 
transforms according to 

whereas other metric independent tensor fields stay in¬ 
variant. It follows that the calculation of the photon 
wave 4-vector described above is gauge invariant. To 
be reassured, let us consider a case study which exem¬ 
plifies that gauge transformations cannot induce photon 
energy shift. Let the “metric perturbation” be gener¬ 
ated by = /(t — x) along with other = 0. Though 
merely gauge, this form resembles a travelling GW in 
the x-direction. Yet, does it influence the energy of co¬ 
moving photons? To check this, consider the resulting 
nonzero metric perturbation components given by 

^XX - 2/*, htx - f- 


(19) 










4 


It follows at once from ®, (0 and noting that hP°‘ = 
—hoa that u! = kx and kx = 0, i.e. no photon energy 
shift, as expected. 

We shall estimate the size of the effect of the pho¬ 
ton energy upshift by GWs from a compact astronomical 
source. To keep our discussion generic it is not necessary 
to analyze any specific theoretical template for gravita¬ 
tional radiation of high post-Newtonian order. Partly 
inspired by the recent suggestion that the core-collapse 
model for GRBs with gravitational radiation may involve 
fragmentation into two in-spiral compact objects we 
proceed with the present simplistic lowest order approach 
to GWs from a binary system. We expect the study of 
such a system to represent similar processes and to pro¬ 
vide good order-of-magnitude guidance on photon energy 
upshift by gravitational radiation from violent astronom¬ 
ical events. 

Consider a binary of black holes of masses mi and 
TO 2 with separation R. Denote the total mass by M := 
mi + m 2 and reduced mass by /i = mim 2 /M. Choosing 
the origin of (a;“) as the centre of mass, the orbit of the 
reduced mass on the x-y plane is given by 0: 

X = R cosilt, y = R sin ilt, z = 0 (20) 

where il = ijM/R?. This yields the nonzero components 
of the quadrupole moment: Qxx = cos^ ilt, Qyy = 
sin^ ilt, and Qxy = fJ,R^ sin fit cos ilt. It then follows 
from m that 

Q = AfiR^n^ sm2ilLo. (21) 

Photons starting from the initial position Lq such that 
Q > 0 above will be subject to an energy upshift in accord 
with the power law II3- The maximum energy upshift 
occurs when Q = Qmax = It is instructive 

to introduce A := m 2 /mi as the ratio of the masses and 
a := R/M as the separation in the unit of the total mass. 
In terms of these dimensionless parameters the maximum 
energy upshift exponent takes the form 

4Aa-5/2 

Qmax = (22) 

It is remarkable that Qmax is independent of the total 
mass of the system. Here the greatest effect is achieved 
with two equal masses, i.e. A = 1. In this case we have 
simply Qmax = Its magnitude falls sharply for 

a > 2 corresponding to a separation larger than the sum 
of the black holes’ radii. Although our calculation above 
is based on the lowest order approach, it has an indicative 
value. In particular, this expression suggests that a linear 
photon energy upshift with distance may be attainable 
for binary black hole mergers with cr = 1. 

The above discussion suggests that the GWs from a 
black hole forming event may interact strongly with pho¬ 
tons from (the vicinity of) the same source, regardless of 


its total mass. For the formation of a solar-sized black 
hole, if we take Q = 1, then, using 113, a photon with 
initial energy 1 eV at Lq = 150 km from the source will 
acquire an energy of up to 1 MeV, i.e. in the gamma 
ray domain, when travelling a distance of only 1 AU. 
This photon may travel a further 15, 000 light years, cor¬ 
responding to a moderate 10% of a galactic diameter, 
to reach the energy of lO^^eV before being scattered by 
GMB. The scattering could initiate a cascade of cosmic 
rays. 

In conclusion, by analyzing the “surfing” of photons 
in GWs emitted by a compact source, we have arrived 
at a universal mechanism whereby photon energy may 
be significantly upshifted. We have demonstrated that 
photons, under some conditions, are energized linearly 
with distance. They are also bunched into structures on 
the GW wavelength scale. This may count for the source 
of energy carried off by photons observed as gamma ray 
events. On applying to other astro-particles such as neu¬ 
trinos, the universality of this mechanism may bear im¬ 
portant implications on further problems in ultra-high 
energy cosmic rays. The progress on research into these 
directions will be reported elsewhere. 
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